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Equations of motion are derived for‘the transient response, to 
shock wave, of a submerged shell with internal structures. A substructuring 
procedure, which does not require calculation of a system stiffness matrix, is 
employed to obtain these equations in a general manner for arbitrary 
internal structures approximated bv finite elements. 
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I INTRODUCTION 

In an earlier paper (Ref. [1]), equations of motion for the transient 
response, to a shock wave, of a submerged shell with internally attached 
structures were obtained by a substructuring procedure. Several quantities 
in these equations, relating to the attachments, were only broadly defined, 
and determined in detail only in an ad hoc manner for specific appendages. 

In this report, matrix expressions for these quantities are defined which 
permit their computation for arbitrary internal structures approximated by 
finite elements. 

The present report addresses the same basic problem as Ref. [1], but 
does not consider the details of the fluid loading on the shell, which may be 
found in Ref. [1], The effects of external loads applied to the internal 
structures, which were not considered in Ref. [1], are included here. 

The substructuring procedure described herein is intended for use in 
transient response problems. Unlike many other substructuring schemes, which 
are used primarily to solve steady-state problems (see, e.g., Ref. [2]), the 
present method does not involve the modes and natural frequencies of the combined 
system (shell and appendages). In addition, a system stiffness matrix is not 
required, since the present technique employs the interaction forces and 
moments at the junctions of the shell and the attached structures. 
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II GENERAL KINEMATICS 


Consider a shell S of arbitrary geometry (see Fig. 1). Let n, 2 and 3 

denote mutually perpendicular directions, respectively normal and tangent to 

the middle surface of the shell. If d , d. and d_ denote corresponding 

n i j 

displacement components of this surface, a 3 by 1 displacement matrix d(s, t) 

may be defined, in which s denotes an arbitrary point on the middle surface, 

t is time, and underlining indicates a matrix; d is assumed positive inward. 

n 

Corresponding to this displacement matrix, the 3 by 1 jth shell mode <£.(s) 
has components (f>^, and <J>^, with which are associated a natural frequency 
w. and a generalized mass 


f T 

U. = md>. d) dA 
J J 
A 

in which m denotes the mass per unit of surface area A of the middle surface 
of the shell and the superscript T denotes a transpose. The shell modes, 
which satisfy the orthogonality condition 


( 1 ) 


j m <t>^ (ftj dA = 0 if k f j (2) 

A 

may be used as expansion functions for the shell displacements, i.e., 

v s 

d(s, t) = l q ,(t) <J> (s) (3) 

k=l b k 


in which the q^^t) are generalized coordinates for the shell and is 
the number of modes of the shell considered. 


Assume that an elastic structure O, entirely within the shell S, is attached 
to S at points Sj(I = 1, 2,..., N^), as shown in Fig. 2. Let 0 be approximated 
by finite elements with interface (attached) nodal points I, coinciding with 


Sj, and interior (unattached) nodal points a(a = 1, 2,...,Ny), corresponding 
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to material points of o not in contact with S. The finite element representation 
of 0 is assumed to include rotational degrees of freedom. 


At each interface point I, certain components of displacement and/or 
rotation of 0 must match corresponding quantities of S, depending on the 
type of support. Thus, at an interface point I, define the vector 


/ = 
%I 




X 

-OI 


(4) 


where the upper partition, an N^. by 1 matrix, contains those physical degrees 
of freedom of 0 which must match corresponding physical degrees of freedom of 

S, and the lower partition, an by 1 matrix, contains the remaining physical 

*) 

degrees of freedom of 0 active at interface point I. 


Since the interior points of o are not connected to S, define the vector 
of physical degrees of freedom of o at an interior point a as 


x u 


X 

-0a 


(5) 


It may be observed in Eqs. (4) and (5), and in what follows, that the symbol - 
above a quantity indicates a quantity associated with physical degrees of 
freedom of 0 constrained to move with S, while the symbol ~ above a quantity 
indicates a quantity associated with physical degrees of freedom of 0 not 


At a simple support, e.g., the displacements of o and S must be equal, 
but the rotations of 0 and S are not equal. Thus, assuming two active 
physical degrees of freedom, the displacement of 0 is located in the 
upper partition of Eq. (4), and the rotation of O is located in the 
lower partition of Eq. (4). 










constrained to move with S. 


If the appropriate definitions of Eqs. (4) and (5) are applied to all 
nodal points of 0, the physical degrees of freedom of a may be arranged in 


the form 


in which the upper partition, an N by 1 matrix, contains the terms from 
all interface points [Eq. (4)] and the lower partition, an N by 1 matrix, 
contains the terms from all interface points [Eq. (4)] and the 
terras from all interior points [Eq. (5)]. 


In the sequel, the elements in the upper partition of Eq. (6) will be 
referred to as "constrained physical degrees of freedom of a", while the 
elements in the lower partition will be referred to as "unconstrained physical 
degrees of freedom of a". The matrix X^ has size N by 1, where 


N = N + N 


is the total number of active physical degrees of freedom of o. 


From Eq. (3), the shell displacements at an interface point I may be 
written as 


- (s r ° E 4 = l q sk (t) 4 (s i ) E Ms 

k=l 


where ^ is a 3 by matrix with elements <{^(Sj) and is a V,, by 1 matrix 
containing the generalized coordinates of the shell. Similarly, the rotations 
of the shell cross section at an interface point 1 may be expressed as 


1 l <sk (t > **<V ; *? Ss 

k = l 









in which ^(Sj) den °tes the 3 by 1 matrix of rotations about the n, 2 and 3 

R R 

axes in the mode ^ and i is a 3 by v matrix with elements (s ) . 

“ i b —k I 

Continuity of appropriate components of displacement and rotation of 
shell and substructure at an interface point I requires that 



( 10 ) 


where D^, an by 6 matrix, is a local rotational coordinate transformation 
which transforms quantities expressed in the S-coordinate system into the 

O-coordinate system. Each row of ^ contains at most three non-zero terms.*> 
Use of Eqs. (8) and (9) permits Eq. (10) to be written in the form 


*01 = k % 


where is an by matrix given by 


C = D - 

-I -I .R 

it 


Application of Eq. (11) to all interface points, with reference to Eq. (6) 
for proper ordering, permits the continuity between shell and substructure 
to be written as 


in which C is an N by v g matrix having the form 
*7 ~ 

The coordinate transformation array in Eq. (10) is used for convenience 
in presentation. Actually, only a 3 by 3 coordinate transformation array 
is required at each interface point, the rows of this array applying to 
both displacements and rotations. 





£= Cl 


Equation (13) expresses the constrained physical degrees of freedom of 0 in 

*) 

terms of the free-free modes of S. 

As suggested in Ref. [3], the response of O may be taken as the super¬ 
position of its static response to the actual motion of its supports and 
the dynamic response with respect to fixed supports, i.e.**\ 


X, = ^ \ 


In Eq. (15), c}> is the N by v fixed modal matrix, where V denotes the 

~o o o 

number of fixed-base modes, and is the by 1 matrix of corresponding 
generalized coordinates. The N by N matrix g is obtained by considering 
the static response of 0 to an arbitrary support motion (support "settlement") 
when the nodal loads corresponding to the unconstrained physical degrees of 
freedom vanish. In such a case, Eq. (15) becomes 


<Vstatic ‘ * <4> 


—0 static 


Equation (13) replaces Eq. (5) of Ref. [1] 

d^Sj , t) = d(s^ , t) 


Equation (15) replaces Eq. (4) of Ref. [1]: 


d^s, t) = l q oi (t) ^.(s) + g(s, t) 
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As shown in Ref. [2], the matrix g (sometimes termed the "constraint 
modes") may be obtained by considering the linear static "force-displacement" 
relation for a: 


K(X ) _ = R 

-0 static — 


in which K is the N by N unconstrained stiffness matrix for the substructure 
and R is the corresponding vector of nodal forces and moments. For the 
determination of £, Eq. (17) may be partitioned to correspond to Eq. (6), 
yielding 


k « 

k „ 

-ii ] 

-12 


k „ 

-21 ! 

—22 



X 


” R ” 


—a 




1 

cH 

i 


0 


l_ _J 

static 

— -* 


since the only nodal loads are the reactions. It then follows from the lower 
partition of Eq. (18) that 


(X ) 

-O static 


‘-22 -21 ( Vstatic 


where the superscript -1 indicates an inverse. A comparison of Eqs. (19) 
and (16) shows that 


8 = -Kol 


-22 -21 


Thus, the elements of £ are the static displacements and rotations (including 
rigid-body motion), corresponding to the unconstrained physical degrees of 
freedom of 0, which result from successive unit displacements and rotations, 
corresponding to the constrained physical degrees of freedom of 0 , at 
attachment points I when all other constrained physical degrees of freedom 
of 0 are set equal to zero. 


The use of Eq. (13) permits Eq. (15) to be written as 


= 4 0 3c + £ C £ g 


(17) 


(18) 


(19) 


( 20 ) 


( 21 ) 
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As suggested in Ref. [2], Eqs. (13) and (21) may be combined into the single 
matrix equation 

X, “ 1 C £ (22) 

• • . . c 

in which X is given by Eq. (6), the transformation matrix T has the form 


1 

In 

i 

ol 

& c i 



and the array of generalized coordinates is given by 



In Eq. (22), is an N by 1 matrix, £ is a V by 1 matrix, with 


c 

and T is an N by V matrix. Equation (22) may be written in subscript 
notation as 

X 0k - j, T {, Oj <k ■ '• 2 . » 

j=l 


(24) 


(25) 


(26) 


. ullV- ' 
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III EQUATIONS OF MOTION OF SHELL 


The kinetic and potential energies of the combination of shell S and 


internal structure 0 may be written as 


1 S ^ n •? 

T= i l y.9c- + T I M. x 2 . 
2 A Vsj 2 £ j Oj 


1 r S 22 1 V V 

v ■ 2 l ml Vj "SJ + 2 V Vi x c 


in which the are given by Eq. (26), the K ^ are the elements of the 
partitioned stiffness matrix of Eqs. (17) and (18), and the M. are the 

elements of the partitioned mass matrix of a from an assumed lumped mass 
*) 

formulation . The mass matrix must be partitioned to correspond to Eq. (6) 


—!1 i ^ 


0 | ™ 2 2 


where the upper left partition, an N by N matrix, and the lower right 
partition, an N by N matrix, are diagonal. 


The equations of motion of the shell may be derived from Lagrange's 
equations: 

idfy + ^ ' 5 n + V « * .V 


If Eq. (26) is substituted into Eqs. (27) and (28), one obtains the 
equations which replace Eqs. (6) and (7) of Ref. [1]: 

^Sj’Vj* = 2 l y i q Si + T 0^oj ,q Sj^ 


V(q Sj ,q aj ) 2 l y i w i q si + V a (q oj’ q Sj ) 


L 
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in which the Q,^. are the generalized forces corresponding to the total, fluid 

pressure p acting on the surface of the shell (see Fig. 1) and the Q are 
1 S J 

the generalized forces from the action of 0 on S, exclusive of the contribution 

derivable from the kinetic and potential energies. As shown in the Appendix, 

the Q are the elements of the matrix 

^ j 

/v t * 

S s = C ! R 

where R, an N by 1 matrix, is the vector of known applied nodal loads acting 
on a (see Fig. 2). 

It follows from Eqs. (26) and (24) that 


3X , 
ok 

3 *0k 


L_t. 

II 

N? 

•• v s 

9 V = 


T. • 
kj 

k = 1, 2,.. 

., N) 


which may be seen to reference terms in the upper and lower left partitions 
of Eq. (23). Thus, the use of Eqs. (27), (28), and (32) leads to 


H n r 

dt ( 3^T7 } = U j q Sj + r kj M k X ak 


N N 

9V 2 v C v 

r- = U.co. q 0 . + ) T, . ) K, . X . 

3 q S.i J 1 Sj k=l kj i=l kl 01 


The substitution of Eqs. (31), (33), and (34) into Eq. (30), with reference 
to Eqs. (18), (23), and (29), results in 


Vi + U J W J q Sj = Q Tj - P j (j = 1. 2 ,..., v s ) 


where the P. are the elements of the vector 
J 


P = P + P 
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whose addends are defined as 


I. C [mjj + li 11 V -12 


Z ■ cV [m 22 ^ + k 21 X, + Z 22 X, - *1 


The bracketed term in Eq. (38) may be recognized as the physical 
equations of motion of a. Hence, the contribution of Eq. (38) vanishes, 
and the equations of motion of S become 


U.n +y # (A) q — Q_. + Q . 
j q Sj j j q Sj y Tj v oj 

in which the generalized forces , corresponding to the forces and moments 

exerted by o on S, are 


Q aj'-j, Vk^alc- j, hi V 


which result from the dynamic reactions 


*ok ■ - <M k x ok + j 2 "ki V 


exerted by a on S 


Using Eqs. (22) to (24), the equations of motion of S may now be written 


in the form 


Equations (40) and (41) are generalizations of specific examples 
given by Eqs. (17), (18); (22), (23); (Al), (A2); and (A12), (A13), 
(A16) of Ref. [1]. 






where 6., is the Kronecker delta and the F,, the elastic forces and moments 
Jk i 

exerted by a on S and vice-versa, are the elements of the N by 1 matrix 


- ■ <*11 + i 12 i) £ % - *12 io So (43) 

The effect of inertia forces appears in the double summation in Eq. (42) 

It may be observed that the equations of motion of the shell have stiffness 
coupling only between shell and substructure. 


In the equations of motion of the shell, the contributions of the inertia 
forces and the elastic reactions of a on S may be given a simple physical 
interpretation by referring to Fig. 3. Figure 3a shows an idealized finite 
element model of a substructure 0 attached to a shell S, and Figs. 3b and 3c 
show the manner in which Eqs. (42) and (43) treat the connection of o to S. 

The effect of inertia forces results from the rigidly attached point mass of 
Fig. 3b, while the effect of the elastic reactions of O on S results from the 
"spring force" of Fig. 3c. 
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IV EQUATIONS OF MOTION OF SUBSTRUCTURE 


The equations of motion of the substructure O may be obtained from 


Lagrange's equations: 


&<,£-> + . ( j .,,2 . v , 

dt dq . dq XJj J ’ O 

in which the are the generalized forces corresponding to the known 
applied nodal loads acting on O. As shown in the Appendix, the . 
are the elements of the matrix 


% = *o * 


It follows from Eqs. (26) and (24) that 


5a,!jA. T c (J>1 ' 2 .v 

3 Vj HjJ k - i +V S » 


which may be seen to reference terms in the upper and lower right partitions 
of Eq. (23). Since the upper right partition of Eq. (23) is a null matrix, 
the use of Eqs. (27), (28), and (46) leads to 

d 3T N „ •• 

dt ( 3q a * = \+N X 0k+N 

9V ” ; » 

3q aj " k=i ak ^ i=i X(7i 

The substitution of Eqs. (45), (47), and (48) into Eq. (44) results in 

k £j_ ^Okj (M k+N X ak+N + 'Sc+N.i X ai “ V = ° 

which may be written in matrix form as 


£ I»i E 22 I m + g "i 2 , i k 22 ] 
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If Eqs. (22) to (24) are substituted into Eq. (50), and the matrix products 
are carried out, there results 

2 ^22 2 + 2 ^22 2 + + —22 &> ± % 


+ ^ B 22 £ C £ S = ^ R 


In view of Eq. (20), which defines the constraint modes, the third 
term of Eq. (51) vanishes. Moreover, as discussed in Ref. [2], the fixed- 
base modes of a are the solutions of the eigenvalue problem 

—22 ( 2>j = Vj ^22 ( Vj 

where 0) is the natural frequency of the jth mode, and hence are orthogonal. 
Thus, Eq. (51) may be simplified to yield the equations of motion of O : 




where and are diagonal matrices whose main diagonal terms are, 
respectively, the generalized masses and generalized stiffnesses of a, and 

<• - £ £ » 22 £ ( 54 ) 

is a by N matrix whose columns are the expansion coefficients of a modal 
series representation of the constraint modes. In subscript notation, Eq. (53) 
becomes ^ - Vg 

V + V Vi ■ ^ Vj \ - G Jk J^klV 1 (i ■ 2 . V (55 > 

It may be observed that the equations of motion of the substructure a have 
inertia coupling only between shell and substructure. 


Equation (55) may be recognized as a generalization of Eq. (A8) of Ref. [1]. 







V 


CONCLUDING REMARKS 


The equations of motion for the transient response of a shell S with 
internal structure a are given by Eqs. (42), (43), (53), and (54). It is 
important to note that no total system (S + o) modal stiffness matrix is 
required in the present formulation, as is required when eigenvalues and 
eigenvectors for the entire system are sought (see, e.g,. Ref. [21) - It 
should also be noted that, in view of the stiffness coupling between S and a 
in Eqs. (42) and (43), the modal mass matrix to be inverted in the solution 
of a transient response problem has size Vg by Vg, determined solely by the 
number of modes of the shell employed in the solution. Formulations with 
inertia coupling between shell and substructure require the size of the modal 
mass matrix to be inverted to increase as the total number of modes V increases. 

The physical stiffness and mass matrices of a, K and M of Eqs. (18) and 
(29), respectively, may be obtained by modifying a finite element computer 
code such as SAPIV (Ref. [4]). The constraint modes may be determined by 
means of Eq. (20), or by utilizing a solution option available with SAPIV. Hence, 
it is possible to evaluate all of the matrices required in Eqs. (42) and (43) 
by using only those rows of K and only those main diagonal terms of M which 
correspond to the constrained physical degrees of freedom of a. Once the 
expansion coefficients of the constraint modes [Eq. (54)] and the modal 
coefficients for the interaction forces and moments [Eq. (43)[ have been 
evaluated, the constraint modes and the matrices K and M may be eliminated from 
the final form of the equations of motion of S and a. This elimination 
and the absence of a system modal stiffness matrix will reduce the 
computer core (central memory) required to solve a given problem. 
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The present formulation will also reduce the amount of computer 
time (central processor time) required to set up the matrices in the equations 
of motion, as compared with that required by a method employing a system modal 
stiffness matrix. By way of example, if the formulation of Ref. [2] is used 
for a transient response problem, the matrices 


t CT 

K T C 

(56) 

CT 

c 


T 

M T 

(57) 


must be evaluated. The computer time required to evaluate Eq. (56) may be 
excessive, especially if the problem under consideration requires many modes 
V and/or many physical degrees of freedom of a, and if K is not well banded. 

It should be pointed out that the rearrangement of physical degrees of freedom 
of O to produce the partitioning of Eq. (18) usually results in a poorly banded 
matrix. For a lumped mass formulation, the computer time required to evaluate 
Eq. (57) is not significant compared with that needed to evaluate Eq. (56). 

Thus, a reduction in the amount of computer time required to set up the 
equations of motion is obtained by the elimination of a system modal 
stiffness matrix. A reduction in solution time is also quite probable. 

The free-free modes of the shell S may be determined by means of a computer 
code such as B0S0R4 (Ref. [5]), while SAP1V may, of course, be used to obtain 
the fixed-base modes of the substructure a . If a consistent mass formulation 
is desired for a, the equations of Sections III and IV may be modified in a 
straightforward manner. However, in such a case, SAPIV, which is based on 
a lumped mass formulation, must be replaced by a computer code employing a 
consistent mass formulation. Although the equations derived in Sections 111 


and IV deal with one piece of internal structure, the present method is easily 
extended to include any number of internal appendages. 
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APPENDIX - GENERALIZED FORCES DUE TO APPLIED LOADS ON SUBSTRUCTURE 

The virtual work due to the known applied nodal loads on substructure o 
may be written as 

6W = (6X^) T R (Al) 

in which 6X is the N by 1 matrix of virtual displacements and R is the 

—a — 

N by 1 matrix of applied loads. In Ea. (Al), both the virtual displacements 
and the applied loads correspond to the unconstiained physical degrees of 
freedom of 0 . From Eq. (21), it follows that 

% -io 6*0 +& £% <A2) 

so that Eq. (Al) may be expressed as 

SW = <5^(^ R) + 5%(C T / R) (A3) 

The virtual work of Eq. (Al) may be written in terms of generalized 
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Since the variations of the generalized coordinates are arbitrary, a 
comparison of Eqs. (A3) and (A6) yields 

~ T -'T ~ 

Q S = £ * R 

aT ^ 

% = - 


(A7) 

(A8) 


The components of generalized force associated with the generalized coordinates 
of the shell S are given by Eq. (A7), and those associated with the 
generalized coordinates of the substructure O are given by Eq. (A8). 
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